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DECLE62.001A S>fcT^ 

A METHOD FOR OBTAINING CONSENSUS CLASSIFICATIONS AND 
IDENTIFICATIONS BY COMBINING DATA FROM DIFFERENT EXPERIMENTS 

RELATED APPLICATIONS 

[0001] This application claims priority of U.S. Provisional application no. 
60/441,910, filed January 22, 2003 and EP application no. EP 03447011.2, filed January 15, 
2003. Each of the above applications are incorporated herein by reference in their entirety. 

BACKGROUND TO THE INVENTION 

[0002] Classification techniques have been used in the field of biology to 
determine relationships and patterns in large sets of data. A means for placing an unidentified 
organism into a classification group which reflects its genotype is important to obtain an 
understanding of the organism and its features, particularly in the treatment and study of 
disease and in the assessment of biodiversity. Information which identifies the closest known 
relation to an unknown disease-causing bacterial or viral strain would enable a physician to 
prescribe the most appropriate course of treatment based on such knowledge. Similarly, 
information which could, for example, enable a group of patients suffering from a disorder 
that has several genetic and phenotypical traits to be categorized into nearest-relation groups 
would enable subpopulations of patients to receive the appropriate treatment and/or enable a 
research scientist to discern the most appropriate set of genes and associated biochemical 
pathways for further study. More particularly, the study of infections and epidemic diseases, 
for example, in hospitals or caused by distributed food or contaminated water, can benefit 
from adequate classification of the causative agents, in a way that the sources of an outbreak 
or infection can be traced and eradicated. The development of DNA chips and microarrays has 
increased the amount of data available for analysis, and increased the need for methods that 
can rapidly and reliably analyze raw experimental data. 
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[0003] Classification techniques have been described in the art. PCT patent 
application numbers WO 01/20536 and WO 01/73602 disclose methods for producing 
hierarchical clusterings from large sets of biological data. PCT patent application WO 
01/45026 discloses a method for displaying data resulting from a consensus classification 
analysis in a visually comprehensible format. 

[0004] While methods to determine classifications have been developed to deal 
with the processing of large amounts of experimental data, the problem remains with the 
accuracy of the classifications. Present methods do not take account of the quality of the 
experimental data from which classifications are made, and hence the classifications so- 
produced are distorted, leading to incorrect analyses which may seriously affect the course of 
a treatment regimen. A further problem with conventional classification methods is that they 
are unable to provide a classification when data is missing, or they ignore the fact that an item 
of data is missing and distort the consensus classification so produced. Distorted and incorrect 
consensus classifications have serious implications for medical research and for the treatment 
of diseases. Thus there are time/cost and health benefits in finding a method which can 
produce accurate consensus classifications. 

AIMS OF THE INVENTION 

[0005] It is an aim of the present invention to provide a method suitable for 
producing accurate classifications and which overcomes the limitations of the prior art. It is 
further an aim of the present invention to provide an apparatus for producing accurate 
classifications and which overcomes the limitations of the prior art. 

SUMMARY OF THE INVENTION 

[0006] One embodiment; of the invention is a method suitable for producing a 
consensus classification of organisms using the data derived from two or more experiments 
performed on said organisms or samples thereof comprising the steps of: 

i) obtaining similarity matrices from the said data, 

ii) producing a composite similarity matrix that is a function of said similarity matrices, and 
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iii) producing a consensus classification from said composite similarity matrix. 

[0007] Another embodiment of the invention is a method as described above 
wherein the function of step ii) comprises averaging the corresponding elements of said 
similarity matrices. 

[0008] Another embodiment of the invention is a method as described above 
wherein each similarity matrix is weighted according the number of experimental characters 
used to calculate said matrix, to arrive at the average. 

[0009] Another embodiment of the invention is a method as described above 
wherein each similarity matrix is weighted by a user defined value to arrive at the average. 

[0010] Another embodiment of the invention is a method as described above, 
wherein said experiments produce product size or retention time results, and wherein the each 
element of each similarity matrix is weighted according to the number of bands or features 
associated with that element, to arrive at the average. 

[0011] Another embodiment of the invention is a method as described above 
wherein said experiments are any of electrophoresis, high performance liquid chromatography, 
gas chromatography, capillary electrophoresis, chromatography, thin-layer chromatography, 
and/or mass spectrometry. 

[0012] Another embodiment of the invention is a method as described above 
wherein the function of step ii) comprises the steps of: 

a) linearizing said similarity data matrices, 

b) averaging the corresponding elements of said linearized similarity matrices of step a) 

[0013] Another embodiment of the invention is a method as described above 
wherein step a) comprises the minimisation of equations: 
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wherein is the number or organisms, samples or genotypes, wherein each technique k results 
in a matrix of pair- wise distance values, so that the . distance value obtained between organism 

d kr 

i and j from technique k is given by d kij , wherein d kij = — — with 
S, =- yy^.. , wherein the consensus distance matrix D f7 is considered as the 

* (p-ixp-2)trft k ' ,J 

unknown true universal distance scale and wherein the goal is to search the consensus 
distances D tj , g k and f k so that d k ij = ^(7^.) and D i} =g k (d k iJ ) hold as true as possible. 

[0014] Another embodiment of the invention is an apparatus suitable for 
performing the methods as described above. 

[0015] Another embodiment of the invention is a computer program comprising a 
computing routine, stored on a computer readable medium suitable for producing a consensus 
classification of organisms using the data derived from two or more experiments performed on 
said organisms or samples thereof according to the methods as described above. 

[0016] Another embodiment of the invention is a device suitable for producing a 
consensus classification of organisms using the data derived from two or more experiments 
performed on said organisms or samples thereof according to the methods as described above. 

DETAILED DESCRIPTION OF THE INVENTION 

[0017] As used herein the term, "classification" refers to the classification of 
organisms, or of the data from experiments performed on said organisms or samples thereof. 
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[0018] The term "consensus classification" refers to a classification of organisms, 
or of the data from experiments performed on said organisms or samples thereof, wherein the 
data is derived from two or more techniques, and wherein said techniques produce one or 
more data types. 

[0019] A consensus classification may show the degree of relationship between 
organisms or samples thereof. Consensus classification may indicate the hierarchical ordering 
within a set of organisms or samples thereof Consensus classifications may provide 
information necessary to construct a dendrogram or grouping by means other than cluster 
analysis. 

[0020] By "organism" herein is meant any animal or plant, any cell, bacterium, 
yeast, phage, virus or prion; it includes organisms of any genus, species, sub-type, biotype, 
phenotype or genotype. 

[0021] By "sample" herein is meant a portion which is intended to represent the 
whole. Non-limiting examples of samples are environmental samples such as soil, water, or 
clinical samples such as blood, sputum, stool, urine. Samples may contain many organisms as 
defined above. Samples may comprise extracts of organisms such as their DNA, proteins, 
glycoproteins, or other quantifiable substances. 

[0022] In one embodiment of the present invention, data derived from experiments 
is classified according to the "data type" that is produced by the result of the experiment. 

[0023] In one embodiment of the present invention, the data type is one or more 
"binary characters' 1 . By "binary characters" is meant the result(s) from an experiment is one of 
two outcomes, said result recordable as a binary character. For example, the outcome of the 
experiment might be either "present" or "absent", "positive" or "negative", "high" or "low", 
"dark" or "bright", recordable as a binary 1 or 0. 
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[0024] In another embodiment of the present invention, the data type is one or 
more "continuous-scale characters 1 '. By "continuous-scale characters " is meant the result(s) 
from an experiment is a value which reflects a magnitude, said result recordable as a 
continuous-scale character such as, for example, a decimal number between -100 and +100. 
This category includes, but is not limited to, data from experiments to measure concentrations, 
kinetic properties, intensities, etc. 

[0025] In another embodiment of the present invention, the data type is one or 
more "multistate or categorical characters". By "multistate or categorical characters" is meant 
the result(s) from an experiment is a state or category, said result recordable as a state or 
category name-tag. This category includes, but is not limited to, for example data from 
experiments which result in a colour (e.g. category names red, blue, yellow, black, green), a 
shape (e.g. category names cube, sphere, cylinder), but also sequence-related genotypic 
properties such as genotypes resulting from Multilocus Sequence Typing (MLST) ,Variable 
Number of Tandem Repeats (VNTR) typing, Microsatelite analysis, DNA conformational 
structure analysis (Single Strand Conformational Polymorphism and heteroduplex 
electrophoresis), Single Nucleotide Polymorphism (SNP) analysis, etc. Binary characters may 
be interchangeable with multistate or categorical characters when the result of an experiment 
is one of two categories. 

[0026] In another embodiment of the present invention, the data type is one ore 
more "product size or retention time characters". By "product size or retention time 
characters" is meant the result(s) from an experiment is the magnitude of one property 
recorded as a function of the magnitude of a second property, recordable as a character. This 
category includes, but is not limited to, data from experiments which result in an 
electrophoresis gel, a high performance liquid chromatography (HPLC) chromatogram, a gas 
chromatogram, a capillary electrophoresis chromatogram, a paper chromatogram, a thin-layer 
chromatogram, mass spectrometry techniques such as MALDI-TOF. 
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[0027] In another embodiment of the present invention, the data type is one or 
more "sequence characters By "sequence characters" is meant the result(s) from an 
experiment is a sequence of characters corresponding to, for example, DNA, RNA and/or 
amino acid sequences, recordable as a sequence of characters. 

[0028] In one embodiment of the present invention, characters of data derived from 
experimental tests, observations or recordings taken from samples of organisms (e.g. intra 
and/or extracellular material from one of more organisms) are collected. 

[0029] In one embodiment of the present invention, said data is recorded as an 
"experimental data matrix". 

[0030] As used herein an "experiment" refers to an experimental action or set of 
experimental actions which leads to an observation or a set of observations that can be 
recorded as a dimensional array of measurements for a single organism, sample or genotype. 
Non-limiting examples of experiments include measuring cell length, determining an 
antibiotic resistance profile, incubating a microtiter plate to measure 96 enzymatic activities, 
running an RFLP electrophoresis pattern, a two-dimensional gel, obtaining a DNA sequence. 

[0031] As used herein an "experimental data matrix" means a set of results from 
one or more individual experiments of the same data type. A non-limiting example of an 
experimental data matrix is given in Figure 1-A, wherein four different experiments which 
produce results of the same data type (experiments 1 to 4) have been performed on samples 
from four different organisms (genotype 1 to 4). While the results are presented as a matrix in 
Figure 1-A, the term "matrix", as used herein, does not limit the data to the matrix format; data 
may be recorded or presented in any format. Non-limiting examples of data formats include 
dimensional arrays, lists, coded data, data. 

[0032] As used herein, the term "similarity value" means a value which indicates 
the relative similarity or distance between two or more samples, organisms, and/or genotypes. 
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The similarity value can be a direct result from the experiment (for example, DNA-DNA 
hybridizations) or can be the result of the comparison of two experimental data sets by means 
of a similarity or distance coefficient, known in the art. Since similarity values and distance 
values can be easily converted into one another, the distinction between both types is only 
morphological. Therefore the term "similarity value" is used throughout this document to 
include both similarity and distance values, except when a distance value is explicitly meant. v 

[0033] As used herein, the term "similarity matrix" means a set of similarity values 
derived from two or more formats including one-dimensional, two-dimensional, and multi- 
dimensional arrays, individual experiments of the same data type. The term "matrix", as used 
herein, does not limit the data to the matrix format; data may be recorded or presented in any 
format. Since similarity matrices and distance matrices can be easily converted into one 
another, the distinction between both types is only morphological. Therefore the term 
"similarity matrix" is used throughout this document to include both similarity and distance 
matrices, except when a distance matrix are explicitly meant. 

[0034] In one embodiment of the present invention, the experimental data matrix is 
used to calculate a similarity or distance between organisms, resulting in a similarity matrix as 
defined above. A non-limiting example of a similarity matrix is provided in Figure 1-B, 
wherein the similarity matrix has been calculated from the experimental data matrix shown in 
Figure 1-A. Similarity matrices may be calculated from experimental data matrices using 
techniques known in the art. 

[0035] It is within the scope of the present invention to use a similarity matrix to 
present classifications between, for example, samples, organisms and/or genotypes. In a non- 
limiting example, a clustering algorithm might be applied to a similarity matrix to produce a 
dendrogram or tree. Figure 1-C provides an example of a dendrogram produced from the 
similarity matrix shown in Figure 1-B. 
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[0036] As used herein, a "composite similarity matrix " refers to a new data matrix 
resulting from the combination of two or more similarity matrices. A composite similarity 
matrix may or may not be derived from a single data type. In one non-limiting example, a 
composite similarity matrix is formed from the combination of two or more similarity matrices 
derived from the same data type. In another non-limiting example, a composite similarity 
matrix is formed by combining similarity matrices derived from different experimental data 
types e.g. a similarity matrix derived from results of a DNA electrophoresis gel of a restriction 
digest with a similarity matrix derived from quantification of mRNA expression levels. 

[0037] One aspect of the present invention is a method that can generate a 
consensus classification from the combination of two or more similarity matrices by 
constructing a composite similarity matrix. 

[0038] The invention encompasses methods which may result in a means to view 
said classifications such as, for example, a dendrogram, a Principal Components Analysis 
(PCA), a Self-Organizing Map (SOM), a Discriminant Analysis (DA), and other known 
grouping techniques. 

[0039] The invention encompasses methods which may result in a means to view 
phylogenetic evolution, mutational evolution or clonal relationships such as, for example, a 
phylogenetic tree, an average distance tree, a minimum spanning tree, or any other graphical 
representation that visualizes clonal or evolutionary relationship. 

[0040] The invention encompasses methods which may evaluate the quality of the 
resulting dendrogram, grouping or display of evolution by means of a general quality score or 
quality score at each branching point (co-phenetic correlation, standard deviations, Jackknife 
and Bootstrap tests) and the invention can use these values to optimize and steer an analysis 
procedure. 
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[0041] The inventors have found that the use of composite similarity matrices to 
produce a consensus classification, leads to a surprising increase in the accuracy of the 
classification so produced. The inventors have further found that the combining of similarity 
matrices according to a method of the invention, including when diverse experimental data- 
types are available, such as, for example similarity matrices derived from DNA 
electrophoresis experiments and similarity matrices derived from enzymatic or metabolic 
activity assay studies, the classification so produced is surprisingly more accurate compared 
with conventional methods of the art. 

[0042] One embodiment of the present invention is a method suitable for 
producing a consensus classification of organisms using the data derived from two or more 
experiments performed on said organisms or samples thereof comprising the steps of: 

i) obtaining similarity matrices from the said data, 

ii) producing a composite similarity matrix that is a function of said similarity matrices, and 

iii) producing a consensus classification from said composite similarity matrix. 

[0043] Another embodiment of the present invention is a method as defined above 
wherein the function of step ii) comprises averaging the corresponding elements of said 
similarity matrices. 

[0044] Another embodiment of the present invention is a method as defined above 
wherein each similarity matrix is weighted according the number of experimental characters 
used to calculate said matrix, to arrive at the average. 

[0045] Another embodiment of the present invention is a method as defined above 
wherein each similarity matrix is weighted by a user defined value to arrive at the average. 

[0046] Another embodiment of the present invention is a method as defined above, 
wherein said experiments produce product size or retention time results, and wherein the each 
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element of each similarity matrix is weighted according to the number of bands or features 
associated with that element, to arrive at the average. 

[0047] Another embodiment of the present invention is a method as defined above 
wherein said experiments are any of electrophoresis, high performance liquid chromatography, 
gas chromatography, capillary electrophoresis, chromatography, thin-layer chromatography, 
and/or mass spectrometry. 

[0048] Another embodiment of the present invention is a method as defined above 
wherein the function of step ii) comprises the steps of: 

a) linearizing said similarity data matrices, 

b) averaging the corresponding elements of said linearized similarity matrices of step a) 

[0049] Another embodiment of the present invention is a method as defined above 
wherein step a) comprises the minimisation of equations: 

/=i j=\ 

wherein p is the number or organisms, samples or genotypes, wherein each technique k results 
in a matrix of pair-wise distance values, so that the distance value obtained between organism 

d k . . 

i and j from technique k is given by d kij , wherein d kij - — — with 
2 p M 

S k = yyrf t;; , wherein the consensus distance matrix ZXis considered as the 

unknown true universal distance scale and wherein the goal is to search the consensus 
distances D ij9 g k and f k so that d k ij = (£>..) and D l} = g k (d k ij ) hold as true as possible. 

[0050] Another embodiment of the present invention is an apparatus suitable for 
performing the methods as defined above. 



12 



[0051] Another embodiment of the present invention is a computer program 
comprising a computing routine, stored on a computer readable medium suitable for producing 
a consensus classification of organisms using the data derived from two or more experiments 
performed on said organisms or samples thereof comprising according to the methods as 
defined above. 

[0052] Another embodiment of the present invention is a device suitable for 
producing a consensus classification of organisms using the data derived from two or more 
experiments performed on said organisms or samples thereof comprising according to the 
methods as defined above. 
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Averaging of similarity matrices 

[0053] In one embodiment of the present invention, a method of obtaining a 
consensus classification comprises the steps of: 

i) calculating two or more similarity matrices from two or more experimental data matrices, 
using methods known in the art, wherein the data type of one experimental data matrix may be 
the same as that of other experimental data matrices in the calculation, or the data type of one 
experimental data matrix may be different from that of other experimental data matrices in the 
calculation, 

ii) calculating a composite similarity matrix by averaging the corresponding elements of the 
respective matrices, i.e. between the same pairs of organisms or samples, and 

iii) calculating a consensus classification from the composite similarity matrix of step ii). 

[0054] According to one embodiment of the invention, a composite similarity 
matrix is calculated by averaging the corresponding elements of the respective similarity 
matrices. This means, for example, in the case of element 1, that sum of element 1 from 
similarity matrix 1, element 1 from similarity matrix 2, element 1 from similarity matrix 3 etc. 
is divided by the number of matrices, so producing an average of element 1 for the similarity 
matrices used in the calculation. 

[0055] In one embodiment of the present invention, a method of obtaining a 
consensus classification comprises the steps of: 

i) calculating two or more similarity matrices from two or more experimental data matrices, 
using methods known in the art, wherein the data type of one experimental data matrix may be 
the same as that of other experimental data matrices in the calculation, or the data type of one 
experimental data matrix may be different from that of other experimental data matrices in the 
calculation, 

ii) calculating a composite similarity matrix by averaging the corresponding elements of the 
respective matrices, i.e. between the same pairs of organisms or samples, wherein each 
element is weighted by user-defined parameters, and 

iii) calculating a consensus classification from the composite similarity matrix of step ii). 
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[0056] In one embodiment of the present invention, a method of obtaining a 
consensus classification comprises the steps of: 

i) calculating two or more similarity matrices from two or more experimental data matrices, 
using methods known in the art, wherein the data type of one experimental data matrix may be 
the same as that of other experimental data matrices in the calculation, or the data type of one 
experimental data matrix may be different from that of other experimental data matrices in the 
calculation, 

ii) calculating a composite similarity matrix by averaging the corresponding elements of the 
respective matrices, wherein each similarity matrix is weighted according the number of 
experimental characters used to calculate said matrix, and 

iii) calculating a consensus classification from the composite similarity matrix of step ii). 

[0057] The inventors found that the averaging and weighting method of step ii) in 
the above embodiments provides a surprising increase in the accuracy of the resulting 
consensus classification. 

[0058] According to one embodiment of the present invention, two or more 
similarity matrices are averaged according to equation [1], wherein Sc is the resulting 
composite similarity matrix, SI is similarity matrix 1, S2 and similarity matrix 2, and Sn is 
similarity matrix n. Equation [1] is an example of combining two or more similarity matrices 
using unweighted averages. 



s x +s 2 + +s n 

Sc = [i] 

n 



[0059] A non-limiting example of the embodiment is provided in Figure 2-A 
wherein similarity matrix SI is derived from 15 experimental character types and similarity 
matrix S2 is derived from 6 experimental character types. The composite similarity matrix so 
formed is calculated according to equation [1]. 



15 



[0060] According to another embodiment of the present invention, two or more 
similarity matrices are averaged according to equation [2], wherein Sc is the resulting 
composite similarity matrix, SI is similarity matrix 1, a is the number of elements used in the 
experimental data matrix to calculate similarity matrix 1; S2 is similarity matrix 2, b is the 
number of elements used in the experimental data matrix to calculate similarity matrix 2; Sn is 
similarity matrix n, p is the number of elements used in the experimental data matrix to 
calculate similarity matrix Sn. Equation [2] is an example of an equation for combining two 
or more similarity matrices using weighted averages. 



c r ^ aS x +bS 2 + + pS n 

6c ~ — ; — : ^ — [2] 

(a + b + + p) 1 J 

[0061] An example of combining similarity matrices using weighted averages 
based on equation [2], is provided in Figure 2-B, wherein similarity matrix SI was calculated 
from an experimental data matrix containing 15 elements and similarity matrix S2 was 
calculated from an experimental data matrix containing 6 elements. Using weighted averaging, 
the composite similarity value Sc is calculated according to equation [2]. 



[0062] The averaging method according to equation [2] provides a surprising 
improvement in the accuracy of the consensus classification so-produced. 

[0063] In another embodiment of the present invention, the values of a, b...p in 
Equation [2] are determined by the user. The averaging method wherein the user determines 
the values of a, b, ...p may also provide an improvement in the accuracy of the so produced 
consensus classification, when conducted by an experienced user who is familiar with the data 
types being analysed. 

[0064] In one embodiment of the present invention, the similarity matrix calculated 
from experimental data comprising characters belonging to the "product size or retention time" 
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data type, for example DNA electrophoresis gels, is calculated by comparing data sets two by 
two. Figure 3-A illustrates a non-limiting example wherein the similarity values are calculated 
from DNA electrophoresis fragment patterns resulting from a restriction digest using a 
restriction enzyme (experiment 1) and using DNA from three different organisms (genotype 1, 
2 and 3). In Figure 3-A, every matching band is considered as one band occurrence, and every 
unmatched band on either pattern is also a band occurrence. For example in comparing 
genotype 1 and genotype 2, genotype 1 contains 3 bands , genotype 2 contains 2 bands, 2 
bands match, the similarity coefficient would be calculated at 2/3 or 66.667%. Thus the 
similarity matrix derived from experimental data of the "product size or retention time" data 
type, is built-up, according to a method of the present invention, by comparing data sets two 
by two. 

[0065] In another embodiment of the present invention, similarity matrices derived 
from experimental data matrices containing characters of the "product size or retention time" 
data-type may be combined to form a composite similarity matrix by taking the number of 
bands or features associated with each element of each similarity matrix for use as a weighting 
constant of that element. 

[0066] A method of the present invention for combining similarity matrices 
derived from experimental data matrices containing characters of the "product size or retention 
time" data type comprises the use of equation [3]. Equation [3] assumes that N product size or 

retention time data types are combined into a composite similarity matrix, wherein 3°* is the 

S e 

value of element e of the composite similarity matrix (Sc), Di is the value of element e of the 

e 

similarity matrix from product size or retention time data type the number of 

S e 

experimental features associated with that element, and D * is the value of element e of the 
similarity matrix from product size or retention time data type i. 
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[0067] Figure 3-B shows a non-limiting example whereby a composite matrix is 
calculated from similarity matrices derived from different DNA restriction digest results. Box 
[1] shows the results of a DNA restriction digest upon the DNA of the organisms of Figure 3- 
A (Gen 1, 2 and 3), using a different restriction enzyme (exp 2). Box [2] shows the similarity 
matrix produced therefrom, using the method described for Figure 3 -A, and Box [2] also 
shows the similarity matrix resulting from Exp 1. The elements or similarity coefficients of the 
composite similarity matrix shown in Box [3] are calculated using equation [3]. For example, 
the composite similarity value between Gen 1 and Gen 2 (83%) is calculated by weighting the 
similarity value between Gen 1 and Gen 2 measured in exp 1 (66%) by the number of bands 
associated with that measurement (3 bands); and weighting the similarity value between Gen 1 
and Gen 2 measured in exp 2 (88%) by the number of bands associated with that measurement 
(9 bands); summing the weighted values (3x66 + 9x88) and dividing by the total number of 
bands (12 bands) to arrive at the composite similarity value (83%). 

[0068] The inventors have found that the combining of similarity matrices 
according to the methods of the invention described above, including when data-types such as 
those from DNA electrophoresis experiments are available, the consensus classification so 
produced is surprisingly more accurate compared with conventional methods of the art. 

[0069] According to another embodiment of the present invention, the averaging of 
similarity matrices to produce a composite similarity matrix as described above may be 
generally applied to similarity matrices derived from any combination of experimental data 
sets or data matrices. It may include the combination of a similarity matrix derived from one 
data type with a similarity matrix derived from another data type. It may also include the 
combination of a similarity matrix derived from one data type with a similarity matrix derived 



18 



from the same data type. It may also include the combination of more than two similarity 
matrices, each derived from a different data type. 

Linearization of similarity matrices. 

[0070] The method described above to create a composite similarity matrix based 
on the weighted or unweighted average of individual similarity or distance matrices works 
well when the expected discriminatory range for both techniques is comparable, and when the 
matrices are complete, i.e. for each experiment there is a similarity value present for each pair 
of entries. However, when two experimental techniques are performed on the same set of 
organisms, and they generate strongly different similarity or distance levels, the composite 
similarity matrix formed using the method disclosed above can be distorted. For example, 
DNA hybridization and 16S rDNA sequencing are techniques whose discriminatory ranges are 
different. Figure 5 A compares the distance matrices from DNA hybridization values and 16S 
rDNA gene sequences. On the scale of 16S rDNA sequence distance (the X-axis), the DNA 
hybridization-based distances occupy a narrow range close to zero distance, whereas on the 
scale of DNA hybridization (the Y-axis), 16S rDNA sequence distances occupy a narrow 
range on the most distant side of the scale. This effect is due to the nonlinear relation between 
both matrices. Both matrices can be considered as zoomed windows on different ranges of a 
hypothetical linear distance scale. Thus the averaging of similarity matrices derived from 
experiments with different discrimination ranges can lead to a distorted composite similarity 
matrix. 

[0071] The averaging of two similarity or distance matrices with different ranges 
wherein one or both matrices contain missing values can even lead to stronger distortion. For 
example, Figure 5B shows two similarity matrices of values from DNA homology and 16S 
rDNA sequence experiments on the same organism (A, B, C and D). As shown, DNA 
homology values range from 88% to 40%, whereas 16S rDNA gene identity ranges between 
98% and 88%. In addition, the DNA hybridization homology matrix has some missing 
elements. In spite of these missing elements, it is clear from the DNA hybridization homology 
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matrix alone that the four genotypes analysed consist of two groups: [A,B] and [C,D]. The 16S 
rDNA sequence identity matrix also suggests the same groupings, although at a different scale. 

[0072] The composite similarity matrix created from these two techniques shows 
averaged values for [AB], [BC], [BD], and [CD] but for [AC] and [AD] it has taken the only 
available values, 90% and 88%, respectively. The resulting matrix provides a distorted view of 
the relationships between these three organisms, as it suggests [AC] and [AD] to be at least as 
closely related as [AB]. The resulting UPGMA dendrogram also depicts a different 
classification as compared to the two dendrograms derived from DNA hybridization and 16S 
rDNA sequence identity individually. 

[0073] One aspect of the invention is a method for constructing a composite 
similarity matrix that combines the information present in the individual matrices in a way that 
the useful information from each of the constituent matrices is optimally preserved, i.e. by 
respecting the particular discriminatory ranges of the techniques applied. 

[0074] In one embodiment of the present invention, similarity matrices are 
combined by averaging to form a composite similarity matrix wherein the constituent 
similarity matrices are linearized prior to averaging. 

[0075] By "linearize" in reference to a similarity matrix herein is meant adjusting 
the value of each element of said matrix such that those values that fall in the window 
considered useful to classification (the range of the technique) are placed on the same linear 
scale as those useful values of one or more other similarity matrices used in the calculation. 

[0076] One advantage of using an averaging method incorporating a linearization 
step for combining similarity matrices is that the similarity range for both experiments does 
not need to be comparable. 
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[0077] Another advantage of an averaging method incorporating a linearization 
step for combining similarity matrices is that the discriminatory depth (taxonomic depth or 
phylogenetic depth) of the methods used to construct the constituent similarity matrices do not 
need to be the same. 

[0078] Another advantage of using an averaging method incorporating a 
linearization step for combining similarity matrices is that the individual similarity matrices do 
not need to be complete. The method allows incomplete similarity matrices, resulting from 
experimental data matrices with missing experiments, to be combined successfully without 
distortion of the composite similarity values towards one of the constituent values. 

[0079] In one embodiment of the invention, a method for linearizing similarity 
matrices according to the invention comprises a mathematical description to the problem of 
linearizing matrices and a solution therefor. One possible mathematical description and 
solution therefor for linearizing similarity matrices according to the invention is disclosed 
below. 

Attributes of the problem 

[0080] In one mathematical description of the problem according to invention, p 
organisms, samples or genotypes, and n experimental techniques applied on these organisms, 
samples or genotypes are considered. Each technique k results in a matrix of pair-wise 
distance values, so that the distance value obtained between organism i and j from technique k 

is given by^ kJJ . In most cases, these distance values are the result of calculating a 
mathematical distance coefficient on the experimental data sets. 
[0081] Non-limiting examples are: 

• Nucleic acid sequences: number of mutations in aligned sequences 

• Banding patterns resulting from electrophoresis of DNA restriction fragments: 
number of different bands, distance-converted Jaccard or Dice coefficients 

• Enzymatic activity profiles: distance-converted product-moment correlation 

• Micro arrays: Euclidean distance. 
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[0082] According to the invention, experimental data matrices obtained from the 
techniques are not necessarily complete; some experiments may not have been performed for 
some organisms or samples, which results in incomplete distance matrices, i.e. distance 
matrices with missing values. Alternatively, experimental data matrices obtained from the 
techniques are complete. 

Mathematical description of the problem 

[0083] The distance scales from the different experiments are normalized by 
dividing them by their corresponding Root-Mean-Square (RMS) values according to 
Equations [4] and [5] below. 

d, . 

where 



Sk= (» wn rX l^j [5] 

(p-l)(/>-2) M JmX 



[0084] The consensus distance matrix D is considered as an unknown "true" 
universal distance scale, with all the individual distance scales for each experiment being 
mapped on that universal distance scale by some nonlinear function according to equation [6], 
wherein each individual experiment k has its own functional dependence (Figure 4). 

d k4 ,*MD 9 ) 9 Vk 9 iJ . [6] 

[0085] In practical cases, this relationship will not be 100% exact (e.g. due to 
scatter on the measurements and because of practical limitations of the experiments). 

[0086] Equivalently, the consensus distances are connected to the individual 
distances by nonlinear function [7], wherein the ideal case - fk 
D ¥ sg k (d kM ) t Vk,i,j [7] 
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[0087] The goal is to search the functions ^* and % k so that these relations hold as 
close as possible. Each function ^* holds information about the range of the experiment k. 

[0088] There are some considerations that put constraints on the functions ^ k as 
follows. 

1. Since identical organisms should have zero distance for every possible technique, 
fk (0) = 0 , and consequently g k (0) = 0 . 

2. On average, the distance d k ij between a pair of organisms should increase if the 

"true" biological distance increases. However, due to experimental errors, statistical 
fluctuation, or imperfections of the technique, this may not be true for particular cases. But, as 
an overall trend, the inventors have found that it should hold for every well-designed 
comparison experiment. Hence, all functions f k should increase monotonically. As a direct 
consequence, g k should increase monotonically as well. 

[0089] Summarizing the mathematical description of the problem, the goal is to 
search the consensus distances ^and the functions Sk and ^*so that - fk(D[j) anc j 
A> = g k (d kiJ ) as as p 0ss jbi e -j^jg can b e p U t i n a m0 re exact, least square sense, by 
minimizing equations [8] and [9]. 

£lteu-MD,)f>V* [8] 



,'=1 >1 
p M 



YL^j-S^)) ,Vk [9] 



(=1 j=\ 



[0090] There are many equivalent solutions to this problem. This reflects the fact 
that only relative differences in distances have a physical interpretation. There is no 
interpretation of the consensus distance values themselves. For example, applying any 
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monotonically increasing function on the consensus distances will hold a new, equivalent 
solution. 

[0091] In the equations, we ignored the fact that some distances kJJ may not be 
known. However, the extension is everywhere straightforward: the summations should be 
modified to exclude the absent values. We kept the simple forms in order not to overload the 
notations and mathematical formulas. 

MATHEMATICAL SOLUTION TO THE PROBLEM 

[0092] In order to further parameterize the problem, one can write the functions ^* 

and ^ k as a linear combination of a number of basis functions according to equations [10] and 

B f 

[1 1], wherein m is the number of basis functions used, ' is the basis functions, and J M and 

are scalar values that are the unknowns of the problem. 

m 

Md) = 1£f k A(d) [io] 

m 

g k (d) = T,g«A(d) [11] 

l=\ 

[0093] One choice for the basis functions (although not the only possibility) are 
power functions according to equation [12]: 

B t {d) = d l . [12] 
[0094] These expansions are subject to four constraints: 

1- /*(0) = o, 

2. £ t (0) = 0, 

3. f k increases monotonically, and 

4. g k increases monotonically. 
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[0095] One can account for 1. and 2. by only using basis functions that fulfill the 
same criterion: " ® . Points 3. and 4. are discussed later below. 

[0096] The problem is now translated to the minimization of equations [13] and 

fa D 

[14], with unknown values M , and ij . 



P i-\ 



\ 2 



1=1 y=i V 1=1 



, V* [13] 
,V* [14] 



[0097] Unfortunately, these equations are not linear in the unknowns: equation 

f D B (D } 

[13] contains terms that are mixed in JkJ and . iJ . Moreover, * iJ may be a nonlinear 

function. As a consequence, one cannot apply the theory of linear least squares optimization to 

this problem. One possible solution would be to apply a general nonlinear iterative 

minimization algorithm (e.g. Levenberg-Marquard). 



[0098] However, in practical cases, it is sufficient to follow a simpler approach, by 

rr D 

first minimizing equation [14] for an optimal solution for 3kJ and iJ , and then minimizing 

7 

equation [13] for a solution for JkJ . The inventors have found that this approach yields good 
results - equation [13] and [14] are to a equivalent to a high degree: they are in fact each 
other's inverse. A solution that is optimal for one equation will be almost optimal for the other 
one. 

[0099] Although equation [14] is linear in the unknown variables, one has to take 
special precautions in order not to arrive at the trivial, perfect but meaningless solution ® iJ ® 
and SkJ . One solution is to follow an iterative approach as follows: 
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1. Set the consensus distances equal to the averaged individual distances according to 
equation [15]. 

^j=\t^ [15] 

2. Minimize equation [14] for the unknown values g k n keeping Z). fixed. This is a 
standard linear least square problem. 

3. . Minimize equation [14] for the unknown values D l} , keeping g k l fixed. 

4. Standardize D y in some way, e.g. by dividing by the total RMS value 

5. Return to 2. until convergence is achieved. 



[0100] Step 2. needs to take into account the fact that Sk must increase 
monotonically. This can be achieved by using a Quadratic Programming technique rather than 

a simple least square fit. Basically, one evaluates the first derivative of & k in a number of 

fixed distance values z = 1 ---^ r j requiring that the first derivative should be non-negative 

everywhere: % k ~ ® , 2 . For example, the points ^ i could be chosen is such a way that 
they are spread in an equidistant way over the whole distance range. A sufficient number of 

points needs to be used to ensure a consistent non-negative derivative. In each point ^ , this 

condition translates into a linear inequality in the unknowns & kJ using equation [16], wherein 

D' B 

1 the first derivative of the basis function ' . 

m _ 

*Zg kJ B;(d,)Z0,Vk,i [16] 

[0101] The minimization of the quadratic form given by equation [14], together 
with these linear inequalities are a standard problem that can be solved by methods known in 
the art, such as Quadratic Programming. 
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[0102] When a solution is obtained for the consensus distances ij \ it is still 

possible to apply any monotonically increasing function c ^ that has c ^ ~" ^ . This does not 
change anything fundamental to the solution and has only a cosmetic interpretation. [JG1] 

FIGURES 

[0103] Figure 1-A, 1-B, 1-C. Example of a conventional cluster analysis based 
upon 4 organisms each characterized by four characters of the same type. 

[0104] Figure 2-A, 2-B, Combining similarity matrices using non-weighted and 
weighted averaging. 

[0105] Figure 3-A, 3-B. Weighted and unweighted averaging of similarity matrices 
obtained from DNA electrophoresis banding patterns. 

[0106] Figure 4. Figure illustrating the differences in active ranges of results of 
three experiments, dl, d2 and d3. 

[0107] Fig. 5 A. Comparison between distance matrices from two techniques with 
different discriminatory ranges. 

[0108] Figure 5B: Example of averaging two similarity matrices with different 
ranges, and containing missing elements. 

[0109] Figure 6. Xbal and Avrll patterns for six E. coli strains belonging to 
different serotypes, according to Example 1. 

[0110] Figure 7. UPGMA dendrogram obtained from Xbal patterns from six E. 
coli strains, according to Example 1 . 

[0111] Figure 8. UPGMA dendrogram obtained from Avrll patterns from six E. 
coli strains, according to Example 1 . 

[0112] Figure 9. UPGMA dendrogram obtained from unweighted average 
similarities of Xbal and Avrll patterns from six E. coli strains. 

[0113] Figure 10. UPGMA dendrogram obtained from weighted average 
similarities of Xbal and Avrll patterns from six E. coli strains. 

[0114] Figure 1 1 . AFLP clustering of Xanthomonas. 

[0115] Figure 12. DNA hybridization clustering of Xanthomonas. 
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[0116] Figure 13. 16S sequence homology clustering of Xanthomonas. 

[0117] Figure 14. Consensus clustering of AFLP, DNA hybridization and 16S 
rRNA sequences of Xanthomonas based upon unweighted average similarity matrix. 

[0118] Figure 15. Consensus clustering of AFLP, DNA hybridization and 16S 
rRNA sequences of Xanthomonas based upon linearized average similarity matrix. 

[0119] Figure 16. UPGMA clustering of two known species and one unknown 
species of genus A. 

[0120] Figure 17. Histon H3 sequence clustering of different members of a 
eukaryote genus. 

[0121] Figure 18. Consensus clustering of fatty acid composition and Histon H3 
sequences. 



EXAMPLES 
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Example 1: Pulsed Field Gel Electrophoresis (PFGE) using different restriction enzymes to 
determine pathotype of E. coli 0157:H7 strains. 

[0122] Six pathogenic strains of E. coli 0157:H7 are analyzed by means of PFGE 
using two different restriction enzymes: Xbal and AvrIL By means of serological tests, the 
strains have been assigned to specific serotypes. The results of the experiments are shown in 
Figure 6 

[0123] The similarities between the strains are calculated using the Dice 
coefficient: 

where N A is the total number of bands in pattern A, Nb the total number of bands in pattern B, 
and N[ab] the number of common bands between patterns A and B. 



A. FFGE-Xbal 

[0124] Using the Dice coefficient on the experimental data for PFGE-Xbal 
depicted in Figure 1, produces the similarity matrix in Table 1. 



STRl 


100.00 












STR2 


66.67 


100.00 










STR3 


72.73 


60.01 


100.00 








STR4 


72.73 


60.01 


100.00 


100.00 






STR5 


57.15 


33.33 


50.00 


50.00 


100.00 




STR6 


0.00 


33.33 


25.00 


25.00 


0.00 


100.00 
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Table '1. Similarity matrix of Xba\ patterns obtained from six E. coli strains. 

[0125] Cluster analysis using the unweighted pair group method with arithmetic 
averages (UPGMA) results in the dendrogram shown in Figure 7. PFGE-Xbal patterns are 
able to distinguish Al serotype strains from the others, as they cluster together. A2 serotype 
strains, however, could not be clustered together as shown by the dendrogram in Figure 7. 
Furthermore, no distinction is possible between strains 3 and 4 from serotype Al. 



B. PFGE-v4vrII 

[0126] Using the Dice coefficient on the experimental data for PFGE-Avrll 
depicted in Figure 6, produces the similarity matrix in Table 2. 



STR 
1 


100.00 












STR 
2 


75.01 


100.00 










STR 
3 


85.72 


57.15 


100.00 








STR 
4 


50.00 


75.01 


28.57 


100.00 






STR 
5 


66.67 


66.67 


50.00 


66.67 


100.00 




STR 
6 


66.67 


66.67 


50.00 


66.67 


100.00 


100.00 



Table 2. Similarity matrix of Avrll patterns obtained from six E, coli strains. 



[0127] Cluster analysis using the unweighted pair group method with arithmetic 
averages (UPGMA) results in the dendrogram shown in Figure 8. PFGE-Avrll patterns are 
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able to distinguish A2 serotype strains from the others, as they cluster together according to 
the dendrogram in Figure 8. Al serotype strains, however, could not be clustered together. 
Furthermore, no distinction is possible between the two strains from serotype A2. 

C. Combined clustering of Xbal - Avrll by unweighted averaging 

[0128] In a third analysis, the two similarity matrices (Xbal and Avrll) are 
combined to create a new consensus matrix by averaging the corresponding values of the two 
matrices using unweighted arithmetic averages. This results in the following composite matrix 
(Table 3): 



STR 
1 


100.00 












STR 
2 


70.84 


100.00 










STR 
3 


79.23 


58.58 


100.00 








STR 
4 


61.37 


67.51 


64.29 


100.00 






STR 
5 


61.91 


50.00 


50.00 


58.33 


100.00 




STR 
6 


33.33 


50.00 


37.50 


45.83 


50.00 


100.00 



Table 3. Composite similarity matrix obtained by unweighted averaging of Xbal and Avrll 
similarities from six E. coli strains. 

[0129] When a UPGMA dendrogram is calculated from this composite matrix, the 
consensus clustering obtained is shown in Figure 9. The linear combination of both techniques 
is able to distinguish every strain from every other. However, serotype Al strains, and in 
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particular serotype A2 strains are still not grouped together. When looking more closely at the 
similarities between the A2 group strains STR5 and STR6, we find the two bands to be 
different in Xbal analysis, whereas all 5 bands are the same in Avrll analysis. This results in 
0% and 100% similarity in the respective matrices (see above). The unweighted average 
matrix consequently shows 50% similarity between the two strains, and the UPGMA 
algorithm does not cluster them together. 



D. Combined clustering of Xbal - Avrll by weighted averaging 

[0130] In a fourth analysis, the two similarity matrices (Xbal and Avrll) are 
combined to create a new composite matrix by averaging the corresponding values of the two 
matrices using weighted arithmetic averages. The weights are determined thus that they 
compensate for the amount of information produced by each pair of experiments. The formula 
used for the consensus similarity SC is given the formula: 
N A S A +N R S* 



S = 



' B B 



N A +N B 



[4] 



Na, Nb being the number of characters in experiments A and B, respectively, and Sa and Sb 
being the similarity in experiments A and B, respectively. Taking weighted averages results in 
the composite matrix shown in Table 4: 



STRl 


100.00 












STR2 


70.59 


100.00 










STR3 


77.79 


58.83 


100.00 








STR4 


63.17 


66.67 


73.69 


100.00 






STR5 


62.51 


53.33 


50.00 


58.83 


100.00 




STR6 


37.50 


53.33 


37.50 


47.06 


71.43 


100.00 



Table 4. Composite similarity matrix obtained by weighted averaging of Xbal and ^Ivrll 
similarities from six E. coli strains. 



32 

[0131] When a UPGMA dendrogram is calculated from this composite matrix, the 

consensus clustering obtained is that shown in Figure 10. Again, every strain is separated from 

every other. In addition, the strains are grouped according to serotype. When looking at the 

serotype A2 strains STR5 and STR6, they share a similarity of 71.43%, which is the weighted 

average based upon 2 bands in Xbal and 5 bands in Avrll: 

(2x0%) + (5xl00%) =?l 
2 + 5 



[0132] PFGE with Xbal only produces two bands per strain, which is far below the 
minimum to obtain a reasonably significant measure of similarity. Five bands as produced by 
Avrll is still a low number, although more reliable than two. Provided that the experiments are 
merged with equal weight per observation (as obtained in formula [4]), a new composite 
experiment is created containing 7 observations, which results in the more reliable clustering 
as obtained. 
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Example 2: Classification of the genus Xanthomonas based upon 3 techniques: AFLP, DNA 
hybridization, and 16S rRNA gene sequencing. 

[0133] A total number of 29 strains belonging to the genus Xanthomonas together 
with the type strain of its closest neighbor genus, Stenotrophomonas, have been analysed 
using three genomic techniques: (i) AFLP (Amplified Fragment Length Polymorphism), an 
electrophoresis technique in which some 30-80 bands are selectively amplified and 
electrophorized after total genomic restriction analysis; (ii) DNA hybridization, a technique in 
which the renaturation rate of a mixture of equimolar amounts of total genomic DNA from 
two organism is measured, to determine the overall homology between the genomes; (iii) 16S 
rRNA sequence analysis, in which a complete or a partial sequence of the 16S ribosomal RNA 
gene is determined, and sequences from different organisms compared by homology. 

A. AFLP 

[0134] In Xanthomonas, as in many other bacterial genera, AFLP is able to 
distinguish between individual strains, and is suitable for strain typing such as done in 
epidemiological and strain authentification studies. Thanks to the large number of bands 
revealed by the technique, the range of discrimination extends to the variety, subspecies or 
even species level. In the present example of Xanthomonas, strains from the same species are 
usually grouped together (except the more heterogeneous species X. axonopodis), whereas 
within the species, the technique is able to distinguish between pathovars (see Figure 11). This 
is reflected in the cases where strains that belong to the same species but constitute different 
pathovars (X. oryzae, X. axonopodis, X. hortorum, X. translucens): these strains have 
consequently lower similarities with each other than strains that are not classified in different 
pathovars (X. fragariae, X. saccharri, X. vesicatoria, X. melonis, X. cassavae, X. codiaei, X. 
vasicola pv. holcicola, X. populi, X. cucurbitae, and X. hyacinthi). 

[0135] However, deeper phylogenies such as the relationships between the species, 
and in particular, the relationship between its adjacent genus Stenotrophomonas, are not 
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reflected, and the technique would falsely suggest that Stenotrophomonas maltophilia is a 
member of the genus Xanthomonas. 

[0136] The obtained similarity range in the AFLP study is between 10.3% and 

96.3%. 

B. DNA hybridization 

[0137] DNA hybridization is a technique which typically differentiates at the 
species level. Different species can be distinguished from each other by low renaturation rates 
whereas within the same species, renaturation rates are usually high. Figure 12 shows a 
clustering of DNA renaturation rates between the same strains as used in the AFLP 
comparison. Unlike AFLP, the technique is not well suited to distinguish pathovars from one 
another: most within-species linkage levels are below the experimental error of the technique, 
which is known to be about 6%. Deeper phylogenies, however, are better revealed than by 
AFLP. This is reflected by the fact that all species, including X. axonopodis, can be separated 
from each other. Furthermore, a deeper group is formed by the species X. translucens, X. 
hyacinthi, and X. sacchari, whereas Stenotrophomonas maltophilia is the most distantly 
clustered strain. 

[0138] The similarity range of the present DNA hybridization study is between 
11.3 and 100%. 

C. 16S rRNA sequence comparison 

[0139] As opposed to most other classification techniques, 16S rRNA sequencing 
is known to reveal deep phylogenetic relationships (genus and below). Figure 13 displays a 
cluster analysis of aligned 16S rRNA sequences from the same strains as described before. 
The dendrogram reveals a strongly different taxonomic structure as compared to AFLP and 
DNA hybridization. Within Xanthomonas, three phylogenetic groups can be found, 
represented by a main core consisting of most species (A), a separate group formed by the 
species X. hyacinthi and X. translucens, affecting monocotyls (B), and a third group formed 
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by X. sacchari (from sugarcane) (C). The latter two groups were also suggested by DNA 
hybridization (Figure 12), but are much more prevalent in the 16S rRNA sequence clustering. 
Distantly separated from Xanthomonas is the Stenotrophomonas maltophilia type strain. On 
the other hand, within-and between-species relationships are not dissolved by 16S rRNA 
sequencing. 

[0140] The similarity range of the present 16S rRNA sequence analysis study is 
between 94.4% and 99.9%. 

[0141] It should be stressed that there is no conflict between the clustering revealed 
by AFLP and DNA hybridization on the one hand, and 16S rRNA sequences on the other 
hand. What actually happens is that each technique has its own window in the phylogenetic 
space, and grouping analyses such as dendrograms should be looked at within that space. In a 
technique such as AFLP, grouping levels below 40-50% should not be looked at, because at 
lower similarity levels, the number of incidentally matching bands becomes statistically too 
important compared to the number of real matching bands, i.e., which are identical. Likewise, 
DNA hybridization does not provide reliable similarity values below 40% because of the high 
experimental error of the technique with too different genomic DNA samples. On the other 
hand, 16S rRNA sequences do not allow taxonomic distances to be calculated from one or a 
few different bases out of 1500. Only when a significant number of bases are different (e.g. 10 
or more), the calculated distances become statistically significant. 

[0142] The obvious problem in studies like this is how to obtain one consensus 
clustering, thereby respecting the taxonomic level of each technique, and thus preserving all 
the information offered by the techniques used. When the similarities of the constitutive 
matrices are averaged to obtain an average similarity matrix, the small but very notable 
differences obtained by 16S rRNA sequencing, are mostly masked by the large similarity 
fluctuations obtained by DNA hybridization and, in particular, AFLP. Indeed, the 16S rRNA 
sequence similarity ranges between 94.7% and 99.9%, meaning a total span of less than 6% 
between the most distant strains. This is much less than the span of the other techniques 
(approximately 90%), and is even less than the experimental error of DNA hybridization. 
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[0143] The logical result is a tree that reflects the species and pathovar 
classification, but does not well reflect the deeper phylogenetic structure of the genus as 
revealed by 16S rRNA sequences (Figure 14). 

[0144] An alternative method, using weighted averages, is not applicable to the 
combination of these techniques, since AFLP patterns in the present study are compared as 
densitometric curves of 2000 values by means of the Pearson product-moment correlation 
coefficient, whereas 16S rRNA sequences are composed of 1500 bases on average. These two 
measures are clearly uncomparable in terms of deriving weighted averages. In addition, DNA 
hybridization values are undefined in terms of numbers of characters. 

[0145] A method according to the invention of linearizing similarity matrices 
which have been derived from different experimental techniques i.e. different data types, has 
been applied to the current data set. The composite similarity matrix resulting from combining 
said linearized similarity matrices according to the invention produces the dendrogram shown 
in Figure 15. 

[0146] The method successfully combines the information delivered by the 
described techniques into a consensus clustering which reflects both the deeper phylogenetic 
relationships of the genera and the species, pathovar, and even strain subdivision of the genus 
Xanthomonas. This study shows that the invention presents an invaluable tool for 
investigating the taxonomic and phylogenetic structure of bacterial taxa, and by extension, of 
all living organisms that are analyzed by a combination of techniques that reveal information 
at different taxonomic depth. Consequently, the method can also be of great use for 
identification purposes, by placing unknown organisms into known classification schemes 
thereby respecting the level of information offered by each technique used. 

Example 3: Grouping members of a eukaryote genus based on fatty acid cell wall 
composition analysis using HPLC and Histone H3 sequencing. 
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[0147] A total number of 15 individuals belonging to a eukaryote genus (taxon not 
further specified) are analyzed by their phenotype using HPLC analysis of cell wall fatty acid 
composition, and by genotype using histon H3 sequencing. The members are sampled from 
different sources: neutral, alkaline, and acid. 

A. Fatty acid cell wall composition 

[0148] HPLC based Fatty Acid Methyl Ester (FAME) analysis is a very sensitive 
technique which usually allows individual organisms to be separated from one another. The 
obtained fatty acid metyl ester profiles are very sensitive to environmental influences such as 
substrate, acidity, temperature. The fatty acid profiles are shown as percentages of total fatty 
acid amount in Table 5. 
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Table 5. HPLC Fatty acid profiles obtained from 15 members of a eukaryote genus 
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[0149] When a UPGMA dendrogram is generated based on Euclidean distance 
calculation between the fatty acid profiles, the grouping obtained is as shown in Figure 16. 
Members of the same species from different sources (neutral, alkaline and acid) are well 
separated from each other, which suggests that the technique allows the origin of a sample to 
be identified based upon its fatty acid profile. Between the species, however, the separation is 
not very clear, especially since A. sp. and A. aberrans do not occur in separate clusters. 

B. HISTON H3 SEQUENCE CLUSTERING 

[0150] Histon H3 sequences of approximately 520 base pairs long were aligned 
and clustered to reveal the deeper relationships between the species. The result is a 
dendrogram depicted in Figure 17. 

[0151] The separation of the two known species and the unknown species is very 
pronounced, and apparently, within A. sylvestris, there is a phylogenetically aberrant member, 
which could not be discovered as such using fatty acid composition analysis. On the other 
hand, it is not possible, using Histon H3 sequence analysis, to discriminate between members 
sampled from sources with different pH. This is not surprising, as in the short term, 
environmental factors will influence the phenotype of the organisms rather than the genotype. 

C. CONSENSUS CLASSIFICATION OF FATTY ACID COMPOSITION AND HISTON H3 
SEQUENCES 

[0152] When a consensus matrix is calculated from the individual similarity 
matrices calculated from fatty acid composition and Histon H3 sequences, the dendrogram 
obtained is as shown in Figure 18. 

[0153] The species subdivision as suggested by the Histon H3 sequences is 
preserved in the clustering, while the phenotypic information as obtained from fatty acid 
composition analysis is also reflected in the tree. Interestingly, the phylogenetically aberrant 
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A. sylvestris member 52441 is not clustered along with the other A. sylvestris members from 
neutral source, but is placed separately based on its sequence divergence. 



